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Abstract. The surface critical behaviour of the semi— infinite one-dimensional 

^ ' quantum Ising model in a transverse field is studied in the presence of an aperiodic 

1} ' surface extended modulation. The perturbed couplings are distributed according to 

[] , a generalized Fredholm sequence, leading to a marginal perturbation and varying 

surface exponents. The surface magnetic exponents are calculated exactly whereas 
C^J ' the expression of the surface energy density exponent is conjectured from a finite- 

size scaling study. The system displays surface order at the bulk critical point, above 
a critical value of the modulation amplitude. It may be considered as a discrete 
Cn i realization of the Hilhorst— van Leeuwen model. 
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1. Introduction 

The influence of bulk quasiperiodic or aperiodic perturbations on the critical properties 
at second order phase transitions has been an active field of research during the last 
years. Different problems were studied numerically on the two-dimensional Penrose 
lattice including the Ising model (Godreche et al 1986, Okabe and Niizeki 1988, 
S0rensen et al 1991), percolation (Sakamoto et al 1989, Zhang and De'Bell 1993) 
and the self-avoiding-walk (Langie and Igloi 1992). In all these cases, no change 
in the critical exponents was observed. Universality was also preserved for three- 
dimensional quasiperiodic systems (Okabe and Niizeki 1990). On the contrary, a 
continuously varying roughness exponent was obtained for interface roughening in 
two dimensions with a modulation of the couplings following the Fibonacci sequence 
(Henley and Lipowsky 1987, Garg and Levine 1987). 

The aperiodically layered two-dimensional Ising model has been also extensively 
studied (Igloi 1988, Doria and Satija 1988, Benza 1989, Ceccato 1989a, 1989b, Henkel 
and Patkos 1992, Lin and Tao 1990,1992a,1992b, You et al 1992, Turban and Berche 
1993) following earlier pioneering works on randomly or arbitrarily layered Ising 
systems (McCoy and Wu 1968, Au-Yang and McCoy 1974). The problem was mainly 
treated in the extreme anisotropic limit where the constant intralayer interaction 
K 1 = J 1 /k B T goes to infinity while the modulated interlayer interactions K 2 (k) go 
to zero, keeping fixed the ratio \ k = K 2 (k)/K^ where K^ is related to K x through 

f Permanent address: Department of Theoretical Physics, University of Szeged, H— 6720 Szeged, 

Hungary 

% Unite de Recherche Associee au CNRS No 155 



o 



>< 



cd 



=ond-mat/9411030 



45 



46 D Karevski et al 

duality (Kogut 1979). In this limit, the physics of the system is governed by the 
one-dimensional quantum Ising model (QIM) with the Hamiltonian 

h-jEH + VK+i]- (i-i) 

fe=l 

Some exact results have been obtained with Fibonacci, Thue-Morse and other 
aperiodic modulations for which the critical behaviour is universal. In other cases 
the Onsager logarithmic singularity of the specific heat was found to be washed out 
(Tracy 1988), like in the random McCoy- Wu model. The structure of the critical 
excitation spectrum and the related conformal aspects have been also explored (Igloi 
1988, Grimm and Baake 1994). 

The situation was recently clarified through the introduction of a relevance- 
irrelevance criterion (Luck 1993a, 1993b, Igloi 1993) generalizing to aperiodic systems 
the Harris criterion for random systems (Harris 1974). The cumulated deviation from 
the average coupling A, A(L) = X^fc=i (^fc — ^)) scales with the chain length L as 
8L^ where 5 is the amplitude of the modulation and u> is the wandering exponent of 
the aperiodic sequence related to the leading eigenvalues of its substitution matrix 
(Queffelec 1987, Dumont 1990). Under a change of the length scale by a factor 
b = L/L', the average thermal perturbation A(L)/L is multiplied by b 1 ^, where v 
is the bulk correlation length exponent, and the amplitude transforms as 

5' = b* /v 6, fc = l + i/(w-l). (1.2) 

The relevance of the aperiodic perturbation depends on the sign of the crossover 
exponent $. In the two-dimensional Ising model with v = 1 the modulation is 
irrelevant when u) < 0, i. e. for the much studied Fibonacci and Thue-Morse sequences, 
marginal when u> = and relevant when to > 0. One expects varying exponents In the 
marginal case and a new type of critical behaviour in the relevant case. 

The three types of critical behaviour were indeed obtained in recent exact 
calculations of the surface magnetization of the QIM (Turban, Igloi and Bcrchc 1994, 
Igloi and Turban 1994, Turban, Berche and Berche 1994). 

In the present work, we study surface aperiodic perturbations generated through 
the Fredhom sequence (Dekking et al 1983) and its generalizations. Such sequences 
lead to a vanishing density of defects in the bulk of the system and the bulk critical 
properties are left unchanged. On the other hand, they induce a marginal surface 
extended perturbation so that the surface critical exponents are nonuniversal, varying 
with the modulation amplitude. Such aperiodic perturbations may be considered as 
discrete realizations of the Hilhorst-van Leeuwen model (Hilhorst and van Leeuwen 
1981, see also Igloi, Peschel and Turban 1993, for a recent review). 

The properties of the generalized Frcdholm sequence are studied in section 2. 
The surface magnetization is calculated exactly in section 3 and the surface energy 
exponent is obtained through finite-size scaling in section 4. The results are discussed 
in the last section. 



2. Generalized Fredholm sequence 

We consider a generalized Frcdholm sequence generated through substitution on the 
three letters A, B and C 
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A -► S(A) =ABCC ■■■ C 

B^S(B) =BCC C ■■■ C (21) 

C -» 5(C) = C C C C ■■■ C 

in 

which is the characteristic sequence of the powers of m. With words of length to = 2, 
one recovers the usual Frcdholm sequence (Dekking et al 1983). 

The substitution matrix M, with entries giving the numbers of A, B and C in 
S(A), S(B) and S(C), then reads 

/ 1 0\ 

M = 1 1 1, (2.2) 

\m — 2 m — 1 my 

with eigenvalues f^ = to, f2 2 3 = 1. The wandering exponent is given by 
lnilo 



ln^ 



= (2.3) 



so that, according to equation (1.2), the Frcdholm modulation is a marginal 
perturbation for the QIM. The numbers of letters of each type at the nth step 
in the inflation process are given by the matrix elements of M n . The asymptotic 
letter densities p A , p B and p c are related to the components of the right eigenvector 
associated with the largest eigenvalue f^ with here 

Pa = 0, P B = 0, Pc = 1, (2.4) 

independent of m. 

In the following, the couplings in the Hamiltonian (1.1), which are distributed 
according to the aperiodic sequence, will be written as X k = \r$ k where A is the 
unperturbed (bulk) interaction and r characterizes the modulation. We associate an 
unperturbed coupling A, i. e. f k — 0, to the letters A or C and a perturbed coupling 
Xr, f k = 1, to B. For example, starting on A with m = 2, one obtains the following 
sequences after n iterations: 
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(2.5) 



C 

C B C C C 

A B BCBCCCBCCCCCCC 

0100010000000 

Equation (2.1) leads to the following relations for the / fe s: 

f mp+ 2=0 (P>0), / 2 = 1, (2.6) 

f m p+ q =° fa = 3, 4,..., to). 
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They can be used to deduce similar recursion relations for the number of perturbed 
couplings rij = YjL=i fk m a sequence with length j, which are obtained by splitting 
the sum over k into m sums over p, giving: 



n ,ni+i =n l+1 + 1 {l> 0), n x = 0, 
n ml+q = n i+i + 1 (« = 2,3,...,m). 



(2.7) 



Iterating these relations, one may check that f k = 1 when k = m l + l (I = 0, 1, 2, . . .). 
For a sequence with length L = m l , 

n L = / = ^, (2.8) 

lnm 

and the asymptotic density of defects satisfies 

Poo = lim ^ = p B = 0, (2.9) 

i. e. the generalized Frcdholm modulation introduces an extended surface perturba- 
tion in the system. 



3. Surface magnetization 

The surface magnetization m s follows from the asymptotic behaviour of the surface 
spin-spin correlation function, lim t _ >oo (cr^(0)cr^(t)), which gives the square of this 
quantity on a semi-infinite system. Writing the correlation function in the basis 
which diagonalizes the Hamiltonian in equation (1.1), m s can be expressed as the 
matrix element (a \ a x | 0) , where | 0) is the groundstate and | a) the first excited state 
of Ti (Schultz et al 1964). These two states become degenerate in the ordered phase 
A > A c as a consequence of long-range order. 

The Hamiltonian can be put in diagonal form (Lieb et al 1961) 

using the Jordan-Wigncr transformation (Jordan and Wigner 1928) followed by a 
canonical transformation to the diagonal fermion operators r\ v . The fermion excitation 
spectrum is obtained as the solution of the eigenvalue problem 

e „W(fc) = -^(fc)-A fe 0„(fc+l) 
ej v {k) = -\ k _ 1 i, l ,{k-l)-^ u {k) 

where the 4> v {k) and ip v (k) are the components of two normalized eigenvectors which 
satisfy the boundary conditions <j) v {Q) — ip l/ (0) = Q. 

Rewriting a x in terms of diagonal fermions with | a) = T]\ \ 0) , it can be shown that 
m s is also given by the first component ^(1) of the eigenvector corresponding to the 
smallest excitation. According to the first equation in (3.2), in the ordered phase where 
e x vanishes, other components of the eigenvector follow from the recursion relation 

1 (fc + l) = -AfcVi(*)- (3-3) 
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The normalization of the eigenvector then leads to the surface magnetization 
(Peschel 1984) 

\ -1/2 

i+ELK 2 • (3-4) 

j=i k=i j 

For the aperiodic system, with \ k = \r' k , this leads to: 

oo 

m s = [S(\,r)}- 1/2 , S(\,r) = Y / \- 2 3 r - 2n >, n = 0. (3.5) 

3=0 

The critical coupling A c generally follows from (Pfeuty 1979): 

3 

Urn l/jVlnA fc =0. (3.6) 

7 — >oo ^ — ' 

fc=l 

Here, A c = r~ p =° =1 keeps its unperturbed value, as expected for an extended surface 
perturbation. 

In order to calculate the sum in (3.5) let us rewrite it as 

oo 

S(\,r) = l + A- 2 +T(A,r), T(A, r) = J^ A~ 2j V~ 2 ™ 3 . (3.7) 

i=2 

The second sum can be splittcd into m parts as 

oo oo 
T(X r) — Y^ \-2(« li + 1 ) r -2™ m ( + i _|_ \^ ^-2(ml+2) r -2n ml+2 _| 



/=() 



(3.8) 






-2(ml+m) -2n m i + „ 



and, using (2.7), the following functional equation is obtained: 

T(A, r) = -^-— [A" 2 - A" 2 " 1 + (A 2m - 1) T(X m , r)] . (3.9) 

A — 1 

This can be iterated to give: 

oo (m— l)m — 1 

T(A,r)=^r- 2 (' +1 )A- 2m;+1 £ A 2p . (3.10) 

The critical behaviour can be extracted, applying a finite-size scaling method due 
to Igloi. Assuming that the surface magnetization displays a power-law singularity 
with a critical exponent (3 S , S(X,r) behaves as t~ 2,3s with t = A~ 2 — A~ 2 , near the 
critical point A c = 1. It can be shown (Igloi 1986) that the sum S L of the first L terms 
in a power series expansion in A~ 2 asymptotically scales like L 2 @ s at the critical point. 
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Figure 1. Spontaneous surface magnetization m s of the Ising quantum chain with 
the m = 2 Fredholm aperiodic modulation as a function of the square of the reduced 
coupling A c /A for different values of the coupling ratio r. The surface transition is 
first-order when r>r c = v2- 
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Figure 2. Spontaneous surface magnetization m s of the Ising quantum chain with 
the m = 3 Fredholm aperiodic modulation as a function of the square of the reduced 
coupling A c /A for different values of the coupling ratio r. The surface transition is 
first-order when r>r c = \/3. 



One easily verifies that the term l = n of the sum in (3.10) contains the powers of 
A~ 2 from m n + l to m n+1 . Cutting the sum at n—l, one collects the contribution to 
the first L — m n terms of the series expansion and, using (3.7), 
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Figure 3. Ising Surface magnetic exponents S and f3' s versus reduced coupling ratio 
r/r c for the generalized Frcdholm sequence with m = 2, 3, 4. The exponent (3' B is 
associated with the approach towards the critical surface magnetization when r > r c . 



S L=mn (l,r) = 2 + 



1 n—1 

2 E 



1=0 



1 - (mr- 2 y 



1 — mr 2 



(3.11) 



Two regimes have then to be considered: when r < r c 
(mr~ 2 ) n whereas it is 0(1) when r>r c . It follows that 



0. 



1 



lnr 



2 In m 



r < r c 

r > r„ 



m, 



\f™"> S L=m n behaves as 



(3.12) 



When r < r c the surface transition is second order as shown in figures 1 and 2. 
The exponent (3 S depends on the modulation amplitude, as expected with a marginal 
perturbation. This dependence is shown in figure 3. The exponent goes to zero linearly 
at r c as (r c - r)/(2r c lnr c ). 

The vanishing of j3 s for r > r c signals the occurence of surface order at the bulk 
critical point. Since the surface is one-dimensional in the corresponding layered two- 
dimensional classical system, it cannot stay ordered when the bulk disorders and the 
surface transition is first-order in this regime as shown in figures 1 and 2. The existence 
of surface order at A c and above is linked to the localization of the eigenvector <t> l . It 
remains normalizable even at the critical point when r>r c . 



The value of the critical surface magnetization m s 
limit n^oo in (3.11), which gives: 



follows from (3.5) taking the 



2r 2 



m— 1 



r > r„ 



(3.13) 
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Figure 4. Variation of the critical surface magnetization m SjC with the reduced 
coupling ratio r/r c for the generalized Fredholm sequence with m = 2, 3, 4 from top 
to bottom. 



The behaviour of m s c is shown in figure 4. It vanishes with a square-root singularity 
at r c . 

The approach to the critical surface magnetization above r c involves another 
surface exponent fi' s such that m s — m s c ^t@ s . It may be calculated, using the same 
method as above, by considering the series expansion S' L (l,r) = S L (l,r)— [m s c ]~ 2 , 
with L = m n , which scales like L~P°. One obtains: 






^. = 2^-1. 

In m 



-n0' 3 



(3.14) 



The exponent f3' s , shown in figure 3, vanishes linearly at r c as (r — r c )/(r lnr c ). 

At the critical value of the modulation amplitude, according to (3.11), the finite 
series expansion gives 



fry — \ 

S L=m n(l,r c ) - 2 + — n ~ In L, 



(3.15) 



so that the surface magnetization vanishes there with a logarithmic singularity 

™ s ( r c)~(-^t)- 1/2 . (3.16) 

The variation of rn s (r c ) with A -2 is shown in figures 1 and 2. 
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4. Surface energy 

The critical behaviour of the surface energy can be studied by considering the finite- 
size behaviour of the matrix element e s = (e | a\ | 0) which does not contain any 
regular part and scales like L~ x " s where x e is the dimension of the surface energy 
density. Here, the state | e) is the lowest two-fermion eigenstate ri\ri\ | 0) of the 
Hamiltonian (3.1). Writing a\ in terms of diagonal fcrmions, one obtains 

e s = ^(1)0 2 (1) - ^(1)^(1) = (e 2 - 6^(1)^(1) (4.1) 

where, in the last expression, we used the relation ip v (i) = — e„ </v(l) which follows 
from the second equation in (3.2) for the surface components of the eigenvectors. 

The matrix element e s on finite systems with size L of the form m n was obtained 
through a numerical solution of the eigenvalue problem 

A fe _ 1 ^(fc-l) + (A 2 fe _ 1 + l)^(fc) + A fe ^(fc + l)=e 2 ^(fc), 

(4.2) 
0„(1) + A^ w (2) - e 2 ^(l), 4> V (L + 1) = 0. 

For small chains (L < 2 15 or 3 9 ), the complete excitation matrix was diagonalized 
while for longer chains (up to L — 2 17 or 3 10 ), equation (4.2) was rewritten as a matrix 
recursion relation 

*-&-(*=f* -1F)U%) 

and varying e v , we looked for the zeros of 4> V (L+1), with 0„(O) =0 and <j) v (X) arbitrary 
in the first vector. Once the zeros corresponding to the two lowest eigenvalues are 
found, the corresponding eigenvectors are normalized in order to evaluate (4.1). The 
exponent x e is then deduced from the slope at large L in a log-log plot. The results 
are shown in figures 5 and 6. 

The behaviour of x e below r c can be deduced from finite-size scaling 
considerations. Low energy excitations scale as L^ 1 and due to the factor e v , the 
1. h. s. of the first equation in (3.2) can be neglected. The leading finite-size 
behaviours of 1 (1) and <j> 2 0-) are the same and follow from (3.4) by cutting the sum 
at j = L. The calculation proceeds like for /3 S and gives 4> 1 (l) ~ 02(1) ~ L^ Xm " with 
x m s = /3 S = 1/2 — In r/ In m. As a by-product, we recover v — [3 s /x ms = 1 for the bulk 
correlation length exponent, as expected for a surface perturbation. Collecting these 
results in (4.1), we obtain: 

e.W-L- 1 - 2 "-, x es =2-2^, r<r c . (4.4) 



When r>r c one finds different size-dependences for the first and second excitations 
as well as for the corresponding components of the eigenvector. Due to surface 
ordering, the first excitation vanishes anomalously as £- 21nr / lnm j j. c . quicker than 
the higher ones with the usual L^ 1 behaviour. Equation (3.3) can still be used to 
calculate the size-dependence of ^(1) which, following the same steps as in the last 
section, gives m s (L) — constant term + L~ Xm s where x' m =f3' s . The term on the 1. h. s. 
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Figure 5. Surface energy density exponent as a function of the coupling ratio r 
for the Fredholm sequence with m = 2. The points correspond to finite-size scaling 
results on chains with sizes of the form 2 n uptoL = 2 10 , 2 12 , 2 17 from top to bottom. 
The heavy line gives the conjectured analytical result. 
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Figure 6. Surface energy density exponent as a function of the coupling ratio r for 
the generalized Fredholm sequence with m = 3. The points correspond to finite-size 
scaling results on chains with sizes of the form 3 n up to L = 3 7 , 3 s , 3 9 from top to 
bottom. The heavy line gives the conjectured analytical result. 



in (3.2) can no longer be neglected for 2 (1)- The size-dependence L 1 / 2 lnr / lnm W as 
obtained numerically. Together, this gives: 



e s (L) ~ L 



-l-x' m /2 



lnr 

In to 



r > r . 



(4.5) 
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One may notice that the anomalous behaviour of the first excitation is in agreement 
with a scaling theory for first-order line transitions involving an irrelevant variable 
and leading to a L l ^ 2Xes size-dependence (Igloi and Turban 1993). Numerical results 
for the surface exponent x e in figures 5 and 6, converge smoothly with increasing size 
towards the conjectured analytic expressions. 



5. Discussion 

As already mentioned in the introduction, the aperiodic surface extended perturbation 
treated here, displays some similarities with the Hilhorst-van Leeuwen model (Hilhorst 
and van Leeuwen 1981) in which the coupling, at a distance k from a free surface, takes 
the form 

^ = A(l + ^). (5.1) 

For the Ising model in 1 + 1 dimensions, the marginal case corresponds to a decay 
exponent y = l. The critical surface order appears at a c = 1/2 and the surface critical 
exponents have the following dependence on the perturbation amplitude: 



Ps = X m s = 2 ~ a > X e s = 2 ~ 2a > a ^ a a 

P's = x 'm s = 2a - 1, x Gg — - + a, a > a c . 



(5.2) 



A comparison with previous results in equations (3.12), (3.14), (4.4) and (4.5) shows 
that the exponents of the aperiodic system are recovered through the following 
correspondence: 

(5.3) 



lnm 

In the Hilhorst-van Leeuwen model, the integrated relative perturbation at a distance 
L is given by: 

f^ = ai2l^alnL. (5.4) 

fe=i fe=i 

On the average, the same logarithmic dependence is found with the aperiodic 
perturbation for which, on a chain with length L = m n according to (2.8), the 
corresponding quantity is 

J2^-± = (r-l)n L = r —±lnL. (5.5) 

*-^ A lnm 

fe=i 

This would suggest the identification a = (r — l)/lnm which, actually, is only valid 
up to the first order in r— 1, i. e. for a weak modulation. The fluctuations around 
the average logarithmic behaviour modify this expression for a stronger perturbation. 
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The correct identification can be obtained by considering the finite-size behaviour of 
the critical coupling which follows from equation (3.6) with 

(5.6) 





ln[A c (L)] = 


= - In r — - 


= 


lnr 
In m 


lni 


for the 


aperiodic system, 


whereas 










ln[A c (L)] = 


1 L 

fe=i 


(4) 


~ — c 



- ( 5 - 7 ) 

in the Hilhorst-van Leeuwen model. 

To conclude, one may notice that an experimental realization of the Hilhorst-van 
Leeuwen model with its l/k v decay of the couplings appears difficult while its aperiodic 
counterpart, studied here in the marginal case, could be more easily obtained (at least 
in three dimensions) using appropriate sequences in multi-layer systems. 
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